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■ 

O ' We present the results of our study of X-ray pulsar Cen X-3 using XMM-Newton observations. 

O ' The light curve and the spectrum for this observations were built and Fe triplet within 6.5 — 7 keV 

^ ' region was detected. The geometric model of relativistic accretion disk [1] for iron emission lines Fe I 

O I Ka, Fe XXV and Fe XXVI in 6.4 — 7.0 keV region was applied. The values of disc inclination, inner 

^ ■ and outer radii of the disc and mass of the central compact object (neutron star) were obtained. 

W ■ Intensity variations of these lines during orbital motion were also detected. The largest variation was 

O ■ detected for Fe I Ka line, that agrees with the results of other authors p]. These results conform 

[ the model in which Fe I Ka line forms in hot plasma of accretion disc and highly ionized iron lines 
■ form in outer regions of binary system. Probably the most interesting feature of Cen X-3 spectrum 

^ is Fe XXV triplet which was found by laria et al. [7] from Chandra data analysis. We did not find 

HH ' this triplet in our analysis of XMM-Newton data and explain its absence by the insufficient energy 

. ' resolution of XMM-Newton instruments. 

^ : 

O ■ Introduction 

U: 

^ ■ X-ray pulsars are binary systems consisting of bright massive star and neutron star. The last one 
draws a matter from a massive bright star and forms an accretion disk around itself. The accretion 

■ disk emits X-ray radiation within the range of energies from hundreds eV to tens keV. X-ray pulsar 
^ I Cen X-3 is a pair of blue giant and neutron star at a distance of 8 kpc from the Earth. The mass 
^ ■ of the neutron star is 1.2 Mq and mass of the giant is 20 Mq. This object was discovered in 1967 
i/^ ] with the help of sub-orbital rocket. X-ray pulses with period of 4.84 sec caused by neutron star 

■ spin were detected by UHURU satellite [6j in 1971. The binary system was optically identified 
I in [8]. The binary has eclipses with orbital period of 2.09 days. Out of eclipse phase the iron 

psj emission line in the spectrum of Cen X-3 is detected at 6.5 keV [TOl. This source was observed by 

^ series of different modern high-energy telescopes: gamma-satellite INTEGRAL within the range of 

^ 16 — 60 keV [5], ground-based Cherenkov telescope Mark-6 at energies above 400 GeV [I] and X-ray 

S< space observatories Chandra and XMM-Newton within the range of 0.2 — 12 keV. Our goal was to 

> . process and analyse XMM-Newton data of Cen X-3. 

c3 ; The data 

In our investigations we used the data obtained by XMM-Newton satellite on January 27, 2001 
during 67 kiloseconds observation. These data are public available on HEASARC (High Energ y 



Astrophysics Scientific Archive Resources Center) database web-page j http:/ / heasarc. gsfc.nasa.gov ) . 
XMM-Newton satellite contains three EPIC cameras (MOSl, M0S2 and PNj that allows us to build 
the images, spectra and light curves for X-ray sources. During this observation two cameras, PN 
and M0S2, worked in imaging mode, and MOSl camera worked in the light curve building mode 
(Timing mode). During observations in Timing mode time resolution and sensitivity of the camera 
increase, which is very useful for obtaining spectra and light curves. In this mode it is possible to 
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build only histogram instead of two-dimensional image. For data processing SAS (Science Analysis 
System) shell ver. 1.52.8 was used. 

During this observation, namely in its end, a solar flare appeared, which led to sharp rise of the 
light curve during 900 seconds. This time interval was excluded, and the resulting data contained 
66 kiloseconds observation. For this last time interval we obtained an image using the data from 
PN camera, which has the greatest sensitivity and time resolution compared to other EPIC cameras. 
Another reason for choosing of PN camera is signiflcant pile-up effect in M0S2 camera. This effect 
lies in the following. If two separate photons arrive during the time which is less than time resolution, 
then only one photon is detected by CCD matrix pixel instead of two and its energy is equal to the 
sum of energies of two photons. The influence of pile-up effect on spectrum is the following: firstly, 
continuum level increases at high energies and decreases at low energies; secondly, a part of photons 
is lost at all. 



Analysis of Light Curve and Spectrum 

In paper [101 the formula for calculation of Cen X-3 eclipses moments was derived taking into 
account the changes of system orbital period. Using this formula we determined the orbital phase 
of Cen X-3 system for the beginning and the end of XMM-Newton observation we consider. These 
phase values are -0.03 and 0.33. Phase corresponds to the mid-eclipse and phases ±0.5 correspond 
to the case when pulsar is in front of the star. The light curve of the observation for PN camera is 
presented in Fig. 1. One can see that X-ray flux significantly increases during post-eclipse phase. 

To build the spectrum we select a circle region from the image with radius equal to 48 arcsec, thus 
the circle touches the image border. From the spectrum of this region the background spectrum was 
substracted. The background spectrum was taken from the rest part of the image. It is worth to note 
that we select the maximum possible resolution of the instrument, 15 eV, as the energy resolution of 
the spectrum. The resulting spectrum was handled by XSPEC software package ver. 12.2. 

We excluded — 0.2 keV range, because XMM-Newton is not sensitive in this energy range, 
and obtained the spectrum, which is presented in Fig. 2. This plot was built without taking into 
account any physical models of radiation. We fit the continuum level by the powerlaw model taking 
into account the photoelectric absorption along the line of sight [2j and the iron absorption limit 
at 7.2 keV. Following [7] we used frozen values for spectral power-law index a = 1.2 and hydrogen 
column density Nh = 1.95 • 10^^ cm~^. We obtained the values of normalization factor Npo = 
(2.91 ± 0.06) • 10~^cm~^ s~^keV~^, absorption limit energy Ef = 7.14 ± 0.03 keV and optical depth 
r = 0.22 ±3. In 6 — 7 keV range of the spectrum one can see 3 iron emission lines. The first line of 
triplet is Ka line, which is emitted by neutral and low-ionised iron atoms. The second and the third 
lines are formed by helium-like ion Fe XXV and hydrogen-like Fe XXVI correspondingly. For more 
detailed study of these lines we rebuilt the spectrum for 6 — 7.5 keV region. We fitted emission lines 
by Gaussian and Lorentzian profiles and found that the first one is more suitable. We obtained the 
central energies, widths and intensities for all lines and uncertainties of these parameters. The iron 
emission triplet can be seen in Fig. 3 and parameters of lines are presented in Table 1. In addition to 
Gaussian profile we used "diskline" model describing relativistic accretion disk p[]. This model allows 
us to estimate the angle between accretion disk and line of sight and the inner and outer radii of the 
accretion disk. We found that the sum of diskline and Gaussian profiles is needed to account the 
finite width of Ka line due to thermal broadening. We obtained the following inclination angles of the 
accretion disk: 60° for Ka line and 55° for Fe XXV and Fe XXVI lines. The difference of inclination 
angle for K\a line and other lines is due to the fact, that it is formed in the outer atmosphere of the 
binary system. According to our model the accretion disk has inner radius of 4500 km and outer 
radius (iDy Fe XXV and Fe XXVI lines) of 405000 km. Our results agree within the errors with 
previous observations [9] according to which the inclination angle is equal to 75° ± 12°. The iron 
triplet with diskline approximation is presented in Fig. 4. 
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Comparison with other observations and conclusions 

Beside XMM-Newton Cen X-3 was observed by other satellites in particular Chandra [7] and 
ASCA [3j. These two satellites obtained Cen X-3 spectra in which iron lines in 6 — 7 keV region were 
detected. We compared iron lines parameters from these observations with our results. 

Chandra observed Cen X-3 on December 30, 2000 during of 48 kilosecoonds at orbital phases of 
0.13 — 0.40. Chandra group determined the equivalet width of Ka line and compared it with the 
results of [3] and [H]. We compared our Ka line parameters with Chandra data during post-eclipse 
phase which partially corresponds to XMM-Newton time of observations. The central energies of the 
line are the same for both satellites. Line intensity for XMM-Newton data is less than for Chandra: 
(0.89±0.06) • 10~'^cm~^s~^ and (1.18±0.19) • 10~^cm~^s~^ correspondingly. This difference lies within 
the errors and is negligible. According to line width is equal to 11.5 ±4.5 eV, which is 5 times 

less than for XMM-Newton observation. We fixed line width at 11 eV and obtained increasing 
of the /d.o.f. value from 1.0126 to 1.2463 and corresponding F-statistics F = 66. This means 
that models with 11 keV and 52 keV are similar with probability of 10~^^. Thus we concluded that 
for XMM-Newton observation intrinsic line width is significantly larger and this effect can not be 
explained by the differences in energy resolution of XMM-Newton and Chandra detectors. Chandra 
observed Cen X-3 when it was in outburst state while the time of XMM-Newton observation (4 
week later) corresponds to quiescent state of system. Thus the differences in X-ray emission can 
be explained by the different states of Cen X-3. The values of continuum level and power-law 
normalization factor for XMM-Newton data is 19 times less that also agrees with quiescent state. 

It is clearly seen from the light curve that the source brightness varies significantly with the orbital 
phase. We considered the changes of the intensity of iron lines with phase. This effect was discussed 
at first time in [3]. The authors examined the changes of Cen X-3 iron lines using ASCA X-ray 
satellite. This satellite observed Cen X-3 during June 24 — 25, 1993 at the following orbital phases: 
-0.31/ -0.27 (pre-eclipse), -0.23/ -0.08 (ingress), -0.07/0.12 (eclipse), 0.14/0.19 (egress). XMM- 
Newton observation interval corresponds to a part of ASCA eclipse interval and the whole egress 
interval. For these two phase intervals we built the source spectra from XMM-Newton observation 
and found the parameters of iron lines (see Table 1). We compared parameters of iron lines for XMM- 
Newton and ASCA in eclipse and egress phases. Taking into account the difference of effective areas 
of satellites detectors we found that continuum levels within 6.5 keV region are equal. Intensities of 
all lines from XMM-Newton observation coincides with ASCA data within errors. We confirmed the 
orbital variations of iron lines. Positions of the lines also coincide within errors. When comparing 
the equivalent widths of the lines we revealed the following. During eclipse phase Fe XXVI line is two 
times broader for XMM-Newton than for ASCA, but two other lines coincide within errors. During 
egress phase all lines from XMM-Newton observation have 1.5 times larger width than from ASCA. 
These differences depend on the features of Cen X-3 accretion state at the time of ASCA observation 
(June 24- 25, 1993). 

Therefore the main comparison results are the confirmation of orbital variations, identified by 
the group of Japanese satellite ASCA |3] and quantitative conformity of lines parameters, namely 
positions and intensities. 

In [7] the discovery of Fe XXV triplet instead of single line was reported. The components of 
triplet were interpreted as forbidden transition, intercombination and resonance lines. In XMM- 
Newton spectrum we can see only single Fe XXV line with the parameters given in Table 1. This 
gives us an evidence that Fe XXV triplet could be resolved only at high state of binary system. 
Also we found Fe XXVI line at 6.9 keV, which was observed by ASCA group but is absent in object 
spectrum from Chandra. We found that lines parameters, namely their orbital variations, position 
and intensity, obtained from ASCA and XMM-Newton observations coincide. In [7] the authors 
explained the lack of Fe XXVI line by the stellar wind decreasing. Taking into account high state 
of Cen X-3 during Chandra observation we conclude that stellar wind rate is not the crucial factor 
that affects X-ray luminosity of the system. 
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Table 1: Parameters of Cen X-3 spectrum. Parameters in columns 2 and 3 — 5 were obtained for the whole energy 
interval (0.2 — 12 keV) and 5 — 8 keV interval correspondingly. Normalization parameter is measured in units of 
10~^ cm~^s~^keV~^ for iron emission lines and in 10~^ cm^^s~^keV~-'^ for power-law model for the whole spectrum. 



Parameter 


Full XMM observation 


Chandra egress 


ASCA eclipse 


ASCA egress 


Norm 


29.1 ±5.6 


43.0 ±0.4 


6.5 ±0.3 


34.0 ±7.6 


Line energy {Ka) 


6386 ± 5 


6394 ± 3 


6377 ±15 


6397 ± 7 


Sigma 


36± 15 


52 ±5 


76 ±21 


49± 15 


Norm 


49 ±7 


89 ±6 


11 ±4 


60± 10 


Line energy (Fe XXV) 


6654 ± 5 


6668 ± 3 


6655 ± 5 


6666 ± 6 


Sigma 


46± 14 


46 ±8 


45± 14 


44± 17 


Norm 


53 ±7 


78 ±6 


24 ±3 


65± 10 


Line energy (Fe XXVI) 


6945 ± 7 


6969 ± 3 


6941 ± 5 


6966 ± 7 


Sigma 


43± 19 


53 ±8 


37± 13 


56± 15 


Norm 


37 ±7 


61 ±5 


21 ±3 


52 ±9 


Reduced 


1.0184 


1.0126 


0.9115 


1.0851 
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Figure 1: Light curve of Cen X-3 for XMM-Newton observation. Time counts from January 1, 1998. 
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Figure 3: Iron emission lines fitted by Gaussian model. 
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Figure 4: Iron emission lines fitted by diskline model. 
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